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Abstract. We show that almost all the zeros of any finite linear combination of inde- 
pendent characteristic polynomials of random unitary matrices lie on the unit circle. This 
result is the random matrix analogue of an earlier result by Bombieri and Hejhal on the 
distribution of zeros of linear combinations of L-functions, thus providing further evidence 
for the conjectured links between the value distribution of the characteristic polynomial 
of random unitary matrices and the value distribution of L-functions on the critical line. 



1. Introduction 

Over the past two decades, there have been many new results at the interface of random 
matrix theory and analytic number theory that can be considered as evidence for the zeros 
of the Riemann zeta function being statistically distributed as eigenvalues of large random 
matrices (GUE matrices or Haar distributed unitary matrices); the interested reader can 
refer to [T7], [12] and [19] for a detailed account with many references, and to [H] for 
the function field framework. Since the seminal papers by Keating and Snaith [Ti] , 
it is believed that the characteristic polynomial of random unitary matrices on the unit 
circle models very accurately the value distribution of the Riemann zeta function (or more 
generally L-functions) on the critical line. This analogy was used by Keating and Snaith 
to produce the moments conjecture and since then the characteristic polynomial has been 
the topic of many research papers, and the moments of the characteristic polynomial 
have now been derived with many different methods, e.g. representation theoretic methods 
(see [S1[IB]), super-symmetry method (see [IZ]), analytic methods (Toeplitz determinant 
methods as explained in the lecture by E. Basor in [17|, orthogonal polynomials on the 
unit circle method [15]) or probabilistic methods ([3]), each method bringing a new insight 
to the problem. Many more fine properties of the characteristic polynomial have been 
estabhshed (e.g. large deviations principle in [8], local limit theorems [16], the analogue 
of the moments conjecture for finite field zeta functions [10], etc.). Moreover, thanks to 
this analogy, one has been able to perform calculations in the random matrix world (whose 
analogue in the number theory world seems currently out of reach) to produce conjectures 
for the analogue arithmetic objects (see [2T] for a recent account). 

There are nonetheless certain results that can be proved in both sides, such as Selberg's 
central limit theorem for the Riemann zeta function and the Keating-Snaith central limit 
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theorem for the characteristic polynomial of random unitary matrices (see |I3]). In fact 
Selberg's central limit theorem can be proved more generally for a wide class of L-functions 
(see [20] and [2]). Roughly speaking, an L-function must be defined by a Dirichlet series 
for Re(s) > 1, have an Euler product (with some growth condition on the coefficients of 
this product), an analytic continuation (except for finitely many poles all located on the 
line Re(s) = 1), and must satisfy a functional equation. Such L-functions are expected to 
satisfy the general Riemann hypothesis (GRH), which says that all the non-trivial zeros 
are located on the critical line, the line Re(s) = 1/2. 

Now if one considers a finite number of such L-functions, satisfying the same functional 
equation, then one can wonder if the zeros of a linear combination of these L-functions are 
still on the critical line. The answer is in general that GRH does not hold anymore for such 
a linear combination even though it still has a functional equation (this can be thought 
of coming from the fact that such a linear combination does not have an Euler product 
anymore). But Bombieri and Hejhal proved in [2] that nonetheless 100% of the zeros of 
such linear combinations are still on the critical line (under an extra assumption of "near 
orthogonality" which ensures that the log of the L-functions are statistically asymptotically 
independent). In this paper we will show that a similar result holds for linear combinations 
of independent characteristic polynomials of random unitary matrices. The result on the 
random matrix side is technical and difficult and besides being an extra piece of evidence 
that the characteristic polynomial is a good model for the value distribution of L-functions, 
the result is also remarkable when viewed in the general setting of random polynomials as 
we shall explain it. The main goal of this article is to show that on average, any linear 
combination of characteristic polynomials of independent random unitary matrices has a 
proportion of zeros on the unit circle which tends to 1 when the dimension goes to infinity. 

More precisely, if f/ is a unitary matrix of order ^ 1, let be the characteristic 
polynomial of U, in the following sense: for 2; G C, 

^u{z) = det {In - zU) . 
From the fact that U is unitary, we get the functional equation: 

$f;(2) = (-z)^det(f/)$^(l/^). 

For z on the unit circle, this equation implies that 

^u{z) = R{z)^{-zYdei{U), 

where R{z) is real- valued (with any convention taken for the square root). The fact that 
$(7 has many zeros (in fact, all of them) on the unit circle can be related to the fact that the 
condition needed for $(/ to vanish is in only unidimensional (i.e. R{z) = for a real-valued 
function R). Now, let {Uj)i<^j^n be unitary matrices of order N, and let {bj)i<^j<^n be real 
numbers: we wish to study the number of zeros on the unit circle of the linear combination 

n 
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If we want that F has most of its zeros on the unit circle, it is reasonable to expect that we 
need a "unidimensional condition" for the equation F{z) = if \z\ = 1, i.e. a functional 
equation similar to the equation satisfied by U. This equation obviously exists if all the 
characteristic polynomials satisfy the same functional equation, i.e. the matrices Uj 
have the same determinant. By symmetry of the unitary group, it is natural to assume 
that the unitary matrices have determinant 1. More precisely, the main result of the article 
is the following: 

Theorem 1.1. Let (&j)i^js£n be a family of (deterministic) real numbers, different from 
zero. For N ^ 1, let 

n 

Pn := J]6j$c/^,,, 
i=i 

where {UN,j)i^j!^n is a family of independent matrices following the Haar measure on the 
special unitary group SU{N). Then, the expected proportion of zeros of Fj^ on the unit 
circle tends to 1 when N goes to infinity, i.e. 

E(|{zGU,F^(^) = 0}|) = iV-o(iV), 

where \{z G '^,Fi^{z) = 0}| is the number of z on the unit circle which satisfy Fn{z) = 0. 

The whole paper is devoted to the proof of this result. Before explaining the strategy of 
the proof, we make a few remarks. 

Remark 1.2. Theorem II. II can be stated as 

limE(h{zeU,F^{z) = 0}\] =1. 

Since the random variable — \{zElJ, Fn{z) = 0}| is bounded by 1, in fact the convergence 

holds in all W spaces for p ^ 1. It also holds in probability since convergence in implies 
convergence in probability. 

Remark 1.3. The fact that we impose our matrices to have the same determinant is similar 
to the condition in [2] of the L-functions to have the same functional equation. Moreover, 
in our framework, the analogue of the Riemann hypothesis is automatically satisfied since 
all the zeros of each characteristic polynomial are on the unit circle. 

Remark 1.4. The fact that the proportion of zeros on the unit circle tends to 1 is a re- 
markable fact as a result about random polynomials. Indeed it is well known that the 
characteristic polynomial of a unitary matrix is self-inversive (that is ctAr-A; = exp(z^)afc 
for some ^ G M, if {ak)o<^k^N are the coefficients of the polynomial). As explained in pQ, 
self-inversive random polynomials are of interest in the context of semiclassical approxi- 
mations in quantum mechanics and determining the proportion of zeros on the unit circle 
is there an important problem. Bogomolny, Bohigas and Leboeuf showed that if the first 
half of the coefficients (the second half being then fixed by the self-inverse symmetry) are 
chosen as independent complex Gaussian random variables, then asymptotically a fraction 
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of ^ of the zeros are exactly on the unit circle. Hence we can say that our result is 
not typical of what is expected for classical random polynomials built from independent 
Gaussian random variables. In our framework, we do not even know the distribution of 
the coefficients and we also know that they are in fact not independent. Consequently the 
classical methods which use the independence of the coefficients (or the fact that they are 
Gaussian if one wants to add some dependence) would not work here. Using general results 
on random polynomials whose coefficients are not independent and which do not have the 
same distribution as stated in |9j, one can deduce that the zeros cluster uniformly around 
the unit circle. But showing that they are almost all precisely on the unit circle is a much 
more refined statement. 

We now say a few words about our strategy of proof of Theorem 11.11 In fact we use 
the same general method as in [2], called the "carrier waves" method, but the ingredients 
of our proof are different, in the sense that they are probabilistic: for instance we use the 
coupling method, concentration inequalities and the recent probabilistic representations of 
the characteristic polynomial obtained in [3]. More precisely, for U G U{N) and t G M, we 
denote by Zu(t) the characteristic polynomial of U taken at e~**, i.e. Zu{t) = $[/(e~**). 
Then we make a simple transformation of the linear combination in order that it is real 
valued when restricted as a function on the unit circle: 

n n 

i=i i=i 

Using the fact that Uj G SU{N), one checks that i^e'^'^/^Zu^iO) is real, and that then 
the number of zeros of Fn on the unit circle is bounded from below by the number of sign 
changes, when 6 increases from to + 27r (with 6q to be chosen carefully), of the real 
quantity given by the right-hand side of the equation above. The notion of carrier waves is 
explained in detail in [2], p. 824-827 and we do not explain it again but we would rather 
give a general outline. The main idea is that informally, with "high" probability and for 
"most" of the values of 6, one of the characteristic polynomials Zu. dominates all the others 
(it is the "carrier wave"). More precisely. Lemma [3.81 implies the following: if 6 depends 
only on and tends to zero when N goes to infinity, then there exists, with probability 
1 — o(l), a subset of [6*0, 6*0 + 27r) with Lebesgue measure o(l) such that for any 9 outside 
this set, one can find jo between 1 and such that log \ Zu.^{9) \ — log \ Zu.{9)\ > 6^/[ogN 
for all j 7^ jo- In other words, one of the terms in the sum of the right-hand side of ([1]) 
should dominate all the others. Moreover, Lemma [3.131 informally gives the following: with 
high probability, the order of magnitude of each of the characteristic polynomials does not 
change too quickly, and then, if the interval [^o^ + 27r) is divided into sufficiently many 
equal subintervals, the index of the carrier wave remains the same in a " large" part of each 
sub interval. Now, in an interval for which the carrier wave index jo remains the same, the 
zeros of Zu.^ correspond to sign changes of e^^^/"^ Zu.^{9)^ i.e. the dominant term of ([1]). 

Then, one gets sign changes of i'^ e^^^^"^ F^^e"^^) , and by counting all these sign changes, 
one deduces a lower bound for the number of zeros of F^ on the unit circle. The main issue 
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of the present paper is to make rigorous this informal construction, in such a way that one 
gets a lower bound — o{N). One of the reasons why the proof becomes technical and 
involved is that we have to take into account two different kinds of sets, and show that 
they have almost "full measure": subsets of the interval [6*0, 9o + 27i) and subsets of SU{N). 

More precisely, our proof is structured as follows. We first give two standard results 
(Propositions [2]T]) and l2.2l) . one on the disintegration of the Haar measure on U (N) (indeed, 
most results on random matrices are established for U{N) and we must find a way to go 
from the results for U{N) to those for SU{N)) and the other one which establishes a 
relationship between the number of eigenvalues in a given fixed arc to the variation of the 
imaginary part of the log of the characteristic polynomial. Then we provide some estimates 
on the real and imaginary parts of the log of the characteristic polynomial (Lemmas EH] and 
13. 2p as well as a bound on the concentration of the law of the log-characteristic polynomial 
(Lemma 13.31) . These estimates and some more intermediary one we establish are also 
useful on their own and complete the existing results in the literature on the characteristic 
polynomial. Then we provide bounds on the oscillations of the real and imaginary parts of 
the log of the characteristic polynomial (Lemma 13. 7p . We then introduce our subdivisions 
of the interval [6o,6o + 27r) and the corresponding relevant random sets to implement the 
carrier waves technique. Finally we combine all these estimates together to show that the 
average number of sign changes of ([T]) is at least N (^1 — O ((log A^)^^/^^)) (the exponent 
—1/22 not being playing any major role in our analysis). 



Notation 

We gather here some notation used throughout the paper. 

U{N) stands for the unitary group of order A^, while SU{N) stands for the subgroup 
of elements U{N) whose determinant is equal to 1. fu{N) and ^su{n) will denote the 
probability Haar measure on U{N) and SU{N) respectively. Similarly we denote by Ku(n) 
and Ksu{N) the corresponding expectations. 

We shall denote the Lebesgue measure on R by A. If a > is a constant and if / is an 
interval of length a, then Aq, will denote the normalized measure ^A on the interval /. 

If n is an integer, we note |1, n] the set of integers {1, ■ ■ ■ , n}. 

If £ is a finite set, we note \S\ the number of its elements. 

For n a positive integer, we note lP^(;v) be the n-fold product of the Haar measure on 

SU{N), and E^j^^-jy^ the corresponding expectation. 

If f/ is a unitary matrix, we note for z G C its characteristic polynomial by ^u{z) = 
det {In — zU). For t G M, we denote by Zu{t) the characteristic polynomial of U taken at 
e-'\ i.e. Zu{t) = $c/(e-^*). 

We shall introduce several positive quantities in the sequel: K > 0, M > and 6 > 0. 
The reader should have in mind that these quantities will eventually depend on A^. Unless 
stated otherwise, N ^ 4 and K is an integer such that 2 ^ K ^ N/2, and M = N/K. In 
the end we will use K ~ N/{\ogNf/^'^ and S ~ {log N)-^/^"^ . 
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2. Some general facts 

In this section, we state some general facts in random matrix theory, which will be used 
in the sequel. 

2.1. Disintegration of the Haar measure on unitary matrices. 

Proposition 2.1. Let P[/(Ar) be the Haar measure on U{N), Fsu{n) the Haar measure on 
SU{N), and for 9 & M, let fsuiN),e be the image ofFsu{N) by the application U h- )■ e*^[/ 
from U (N) to U (N) . Then, we have the following equality: 

™ /ON 
^SU{N),e — = ™(7(Ar), (2) 



i.e. for any continuous function F from U{N) to R+, the expectation Ksu{N),eiF') of F 
with respect to 'Fsu{N),e is measurable with respect to 9 and 

f.27r 



Jo 



Proof. One has 

^suiN)AF) = [ F{Xe'')dFsu(N){X), (3) 



which, by dominated convergence, is continuous, and a fortiori measurable with respect to 
9. By integrating ([3]) with respect to 9, one sees that the proposition is equivalent to the 
following: if [/ is a uniform matrix on SU{N), and if Z is independent, uniform on the unit 
circle, then ZU is uniform on U{N). Now, let A be a deterministic matrix in U{N). For 
any dec such that rf"" = det(A), one has Ad e SU{N), and then ZUA = {Z/d){UAd), 
where: 

(1) UAd follows the Haar measure on SU{N) (since this measure is invariant by mul- 
tiplication by Ad e SU{N)). 

(2) Z/d is uniform on the unit circle (since d, as det{A), has modulus 1). 

(3) These two variables, which depend deterministically on the independent variables 
A and Z, are independent. 

Hence ZUA has the same law as ZU, i.e. this law is invariant by right-multiplication by 
any unitary matrix. Hence, ZU follows the Haar measure on U{N). □ 

2.2. Number of eigenvalues in an arc: The result we state here relates the number 
of eigenvalues of a unitary matrix on a given arc to the logarithm of its characteristic 
polynomial. For U G U{N) and t G M, we denote by Zjj{t) the characteristic polynomial 
of U taken at e~**, i.e. Zu{t) = $;7(e~**). Moreover, if e** is not an eigenvalue of U (which 
occurs almost surely under Haar measure on U{N), and also under the Haar measure on 
SU{N), except for e** = 1 and iV = 1), we define the logarithm of Zu{t), as follows: 

N 

logZ^(t):=5^1og(l-e^(^^-*)), (4) 
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where 6i, . . . ,6^ is the sequence of zeros of Zu in [0, 27r), taken with multiphcity (notice that 
the eigenvalues of U are e*^^ , • • • , e*^^ ) , and where the principal branch of the logarithm 
is taken in the right-hand side. We then have the following result, already stated, for 
example, in |H]: 

Proposition 2.2. Let ^ s < t < 2tt, and let us assume that s andt are not zeros of Zu- 
Then, the number of zeros of Zjj in the interval {s,t) is given as follows: 

^ N 1 

he,eis,t)} = ^it-s) + - (ImlogZuit) - lm\ogZu{s)) . (5) 

k=l 

Proof. It is sufficient to check that for all 6 G [0, 27r)\{s, t}, 

7il{eeis,t)} = ^ + Imlog (l - e'^'-'^) - Imlog (l - e'^^-^)) . 

Now, for V G (0, 27r), 

1 - = e™/2(e-W2 _ gW2) _ _2zsin(i;/2) e™/^ = 2sin(t;/2) e'(-— 
Now, sin(t;/2) > and {v - n)/2 G (-7r/2,7r/2) and hence 

Imlog (l-e^") 



V — TC 



since we take the principal branch of the logarithm. Now, for 6 G [0, 27r)\{s, t}, 6 — s + 
27rIL|0<s} and 6 — t + 27rl|0<i} are in (0, 2tt), which implies 



Imlog (1 - e^(^-*)) - Imlog (l - e^^^"")) 



t-n + 27r]I{e<t} 6 - s - it + 27rl{e<s} 



S — t , X 



and then Proposition 12.21 □ 



3. Proof of Theorem 11.11 

3.1. Conventions. All the random matrices we will consider are defined, for some N ^ 1, 
on the measurable space (A^Ar(C), J-"), where J-" denotes the Borel a-algebra of A^Ar(C). 
The canonical matrix, i.e the random variable from (A^iv(C), J-") to A^Ar(C) defined by 
the identity function, is denoted X. Moreover, we denote by ^^uiN) the expectation under 
P[/(Ar), the Haar measure on U{N), and by Ksu{n) the expectation under fsu{N), the Haar 
measure on SU{N). For example, if F is a bounded, Borel function from A^Ar(C) to M, 

EsuiN)[F{X)] = [ F{M)dFsu(N){M). 

JMn{<C) 
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3.2. An estimate in average of the logarithm of the characteristic polynomial. 
Lemma 3.1. There exists a universal constant Ci > such that for all N ^ 2, and A ^ 0, 



Proof. For all A ^ 0, 



2-IT 



SU(N) 



\ogZx{e) ^AVl^j— ^ e-^A^ ¥.suiN){e 



\\\ogZxm 



de 

2^ 



Using the inequality el"'"*"'^! ^ ^a+b ^a-b ^-a+h ^-a-b ^ valid for all a, 6 G M, and writing 
the right-hand side of this inequality as 4E for B and B' being two independent 

Bernoulli random variables independent of U such that P (i? = 1) = 1— P (5 = —1) = 1/2, 
we have: 



27r 



P, 



SU{N) 



(|logZxW| ^Av^b^)^ 



\{B Rc log Zx (0)+B' Im log Zx (0) ) 



We now use the fact ([13] and [7], p. 16) that for Sjt G C such that Re(s + it) and 
Re(s — it) are strictly larger than —1: 



U{N) 



Ke\ogZx{0)+t ImlogZx(O) 



) 



G (1 + ^) G (1 + ^) G (1 + AT) G (1 + AT + s) 



G (1 + AT + iti*) G- (1 + + ^) G (1 + s) 
where G is the Barnes G-function, defined for all 2; G C, by 



(6) 



00 

G{z + 1) := (27r)^/V[(i+^)^'+^]/^ JJ (^1 + 



■z+{z'^/2n) 



n=l 



7 being the Euler constant. The function G also satisfies the functional equation G(2; + l) 
T{z)G{z). 

In other words, one has 



s RelogZx(0)+t ImlogZx(O)^, _ (l + ''t'^) ^ + Ar(g^+*^)/4 

G(l + s) 



E(7(V) (e'' '^<'l°g^^(0)+* ImlogZx(0)j 

where, by the classical estimates of the Barnes function, 

^~{s-+t-)/^ G{l + N)G{l + N + s 



G 



G 



N,s,t 



G{1 + N + ^)G{1 + N+^) 
tends to 1 when A^ goes to infinity, uniformly on s and t if these parameters are bounded. 
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For any sequence {Xn)n^i such that Xn G [0,1/2], one has (taking s = XnB and 
t = XnB'): 



Eu{N) i e^^ ^"^ ^ 



where 



M(A 



N 



: M(Ajv)A^ 
G (1 + A^^) G (1 + A;v^) 



2 



^(l + AArE) 



-G 



N,\i^B,\f^B' 



Since the function G is holomorphic, with no zero on the half-plane {Re > 0}, and since 
Gn,xb,xb' tends to 1 when goes to infinity, uniformly on A G [0, 1/2], the quantity M(A) 
is uniformly bounded by some universal constant c' > 0, for A G [0, 1/2]. Hence, 



E 



U(N) 



^gAAr(B RelogZx(0)+B'ImlogZx(0))j ^ c'iV^, 



for going to infinity, which implies: 



FsuiN) ( \ogZx{e) ^ A^\^) — ^ 4c'e 
Now, taking Ajv = (1/2) A (A/^log N) gives 

^su{N) ( log^x(^^) ^ ^VlogiVj — ^ 4c'e- 



/^-AjvA^/I5ilV+(A2 log N)/2 



< 4^/g-[(v^I^/2)AA](A/2)_ 



□ 



3.3. An estimate on the imaginary part of the log-characteristic polynomial. 

From the previous result, we obtain the following estimate for the imaginary part of the 
log-characteristic polynomial: 

Lemma 3.2. There exists a universal constant > such that for all N ^ 2, 0, and 

e eR, 



P 



SU(N) 



Imlog Zx{9) 



^ AVlogiV ^ c\e 



log AT ^ 



Proof. We use here the probabilistic splitting established in [3] which shows that (see also 
|3] for an infinite- dimensional point of view), for any U eU (N), there exists, for 1 ^ j ^ N, 
Xj on the unit sphere of , uniquely determined, such that 



U = R{xn) 



R{xn-i] 




R{xn-2) 









R{xi] 








(7) 



where R{xj) denotes the unique unitary matrix in U{j) sending the last basis vector Cj of 
C-' to Xj, and such that the image of Ij — R{xj) is the vector space generated by cj — Xj. 



10 



Y. BARHOUMI, CP. HUGHES, J. NAJNUDEL, AND A. NIKEGHBALI 



Moreover, the characteristic polynomial of U{N) is given by 

N 

Zu{0) = l[{l-{xj,e,)), 

i=i 

and then its logarithm is 

TV 

logZ^(0) = ^log(l-(x„e,)), (8) 

i=i 

when 1 is not an eigenvalue of U, taking the principal branch of the logarithm on the right- 
hand side. Notice that the determination of the logarithm given by this formula fits with 
the definition involving the eigenangles (j4]). Indeed, the two formulas depend continuously 
on the matrix U, on the connected set {U G U{N), 1 ^ Spec(?7)}, and their exponentials 
are equal, hence, it is sufficient to check that they coincide for one matrix U, for example 
-In (in this Cj for all j and the two formulas give A^log2). 

If U follows the uniform distribution on U{N), then the vectors {xj)i^j<^N are indepen- 
dent and Xj is uniform on the sphere of C-' . The determinant of U is equal to the product 
of the determinants of R{xj) for 1 ^ j ^ A^, and since R{xi) is the multiplication by xi 
on C, one has 

N 

det(f/) = Xi Y\_^j{xj), 

i=2 

where Tj is a function from to the unit circle U. From this, let us deduce that under 
the measure fsu(N),e- 

(1) The vectors {xj)2^j^N are independent, Xj being uniform on the unit sphere of . 

(2) The value of xi G U is uniquely determined by the determinant det{U) = e*^^, 

N 

x,=e^'^'l[[r,{x,)]-\ 
i=2 

Indeed, let IP5^(^) q be the probability measure on the image of SU (N) by the multiplication 
by e*^, under which the law of {xj)i^j^N is given by the two items above. This probability 
measure can be constructed as the law of the random matrix U given by the formula (^^, 
where {xj)if^j^N are random vectors whose joint distribution is given by the items (1) and 
(2) just above. We now have to prove that fsuiN),9 = ^'su(N)e- ^^^^ notice that 

the joint law of (a;j)2^j<gAr, under the probability measure ^'su{n) (^oes not depend on 9. 
Hence, under the averaged measure 

the vectors {xj)2<^j^N still have the same law, i.e. they are independent and Xj is uniform 
on the unit sphere of 0. Moreover, conditionally on {xj)2^j<^Ni Xi = e"^^^ Y[j=2i^j(-'^j)]~^ ^ 
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where 9 is uniform on [0,27r). Hence, {xj)i^j^N are independent, Xi is uniform on U, and 
then Xj in uniform on the unit sphere of for all j G {1, . . . , A^}, which implies 

2" de ^ /"^"^ de 

^su{N),e ^ - ^u(N) - ^su{N),e 

Now, 'Psu{n),2tt/n is the image of '\Psu(n) by multiplication by e""^^^^ In, which is a ma- 
trix in SU{N): the invariance property defining the Haar measure Fsu{n) implies that 
'^suiN),2Tr/N = ^su{N), and then 9 i— i- ^su(N),e is (27r/iV)-periodic. It is the same for 
9 I— 7- "P'suiN) 6' since the values of xi, . . .xn involved in the definition of ^'su{n) e ^'^ '^'^^ 
change if we add a multiple of l-n jN to 9. Hence, 

2-/^ , Nd9 _ /-^-/^ Nd9 

Now, let F be a continuous, bounded function from U{N) to M. By applying the equality 
above to the function U i— )■ F{U)l^^^^_jj^^^iMe eg/jj) for an interval / C [0,27r/A^), one 
deduces with obvious notation that: 



d9 r„ ,„.d9 
L/T 



^'su{N),9i^^TT\ = / ^su{N),e{F)—, 



where \I\ is the length of I. Now, by definition of fsuiN),e and ^^^(tv) ^^^^ measure 

SU{N) 



is the image of fsuiN) by multiplication by e*^, and the second measure is the image 



/ giTve Q \ 

of ^'su(N) right multiplication by the matrix ( g j 1 • Hence, by continuity 

and boundedness of F, and by dominated convergence, 'Ksu{n),9{F) and K'gjj^j^^ g{F) are 
continuous with respect to 9. By considering a sequence (/r)r>i of intervals containing a 
given value of 9 and whose length tends to zero, one deduces, by letting r — )■ oo, 

^'su{N),ei^) = ^suiN),e{F). 

We now get the equality T'suiN),e = ^'su{n) 6" then the law of {xj)i<^j<^N under fsuiN),e 
described above. 

Hence, the sequence (a;j)2<gj^Ar has the same law under T'su{N),e and P;7(7v)- We now use 
this fact to construct a coupling between these two probability measures on the unitary 
group. 

The general principle of coupling is the following: when we want to show that two 
probability distributions Pi and P2 on a metric space have a similar behavior, a possible 
strategy is to construct a couple {U, U') of random variables defined on the same probability 
space endowed with a probability P, such that the law of U under P is Pi, the law of U' 
under P is P2, and the distance between U and U' is small with high probability. In the 
present situation, we take {x'j)i<^j^N independent, x'^ uniform on the unit sphere of for 
all j G {!,..., N}, and we construct, by using ([7]), a random matrix U' following P(7(Ar). 
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Then, we do the couphng by taking xj := x'j for 2 ^ j ^ N and 

N 

which gives a random matrix U following '^su{N),e- From the fact that Xj = x'j for j ^ 2 
and the equation ([H]), we get the following: 

logZf/(0) - hgZu'iO) = log(l - a;i) - log(l - x[), 

and in particular, 

\lm\ogZu{0) - ImlogZ[//(0)| ^ vr. 

Now, ioT B := lA- 



one gets: 



P 



SU{N) 



Im \ogZx( 



= P 
^ P 
= P 



^ AVlogiv) = P5t/(7V),e (I ImlogZx(0)| ^ AVlogiV 
ImlogZc/(0) ^Av/logiv) 
ImlogZt//(0) ^ A^/hgN -71 
ImlogZx(O) ^^v^logiV 

(|lmlogZx(^)| ^50^)^ 



U{N) 
P 



SU(N) 



/•27r ^fl 

P5C/{JV) ( logZx(^^) ^ 5v/toiiVj — 



^ cie 



Now, if S ^ 



A 



2 V 2 



V\ogN\ / ^2 1 



TT 



22V yiogiV 

€ ^ \ 

2 2 21ok2 



52 



TT 



2 21ogA^ 



V\0gN\ , TT 



+ - 



TT 



2 21og2 



jf ^ > yiogA^ 



2 ' 



2" 



VlogiVX AyiogiV VIogiV 



4 



5 + 



TT 



VIogiV 



5 A/log TT 

4 ^4 



B 
2" 



5 A 
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Hence, we get Lemma [3.2[ with 

/ 2L_| 

C\ = Cie2^21og2 



□ 



3.4. Bound on the concentration of the law of the log-characteristic polynomial. 
Lemma 3.3. For N ^ A, 9 e [0, 2n), xqER and5 E (0, 1/2), one has 
¥su(N)[\hg\Zx{e)\ ~ xo\ ^ ^ C51og(l/5), 

where C > is a universal constant. 
Proof. The proof of Lemma 13.31 needs several steps. 
Sublemma 3.4. For j ^ 1 integer, s,t G M, let us define 

Then, 

(1) For s^ + t^^ 8f, \Q{j, s,t)\^ max (^1 



(2) Forf<:s^ + t^<:8f, \QU,s,t)\^l. 

(3) Fors^ + t^^f, \Q{j,s,t)\ ^ e-(^'+*')/io^''. 



Proof. One has: 



1-^ + tt/j 



If + ^ 8j^, it is immediate that the numerator has a smaller absolute value than the 
denominator, i.e. \Q{j,s,t)\ ^ 1. Moreover, 

and in the case where + ^ j^, one deduces 



\Q{j,S,t)\'^l 



32f 
and then 



14 Y. BARHOUMI, CP. HUGHES, J. NAJNUDEL, AND A. NIKEGHBALI 

Now, if s^+t"^ ^ 8j^, the numerator in ([9]) has a larger absolute value than the denominator, 
and then \Q{j,s,t)\ ^ 1. Moreover, since (s^ + t^)/?>f ^ 1, 



\ 2 



64 ■ 1 + QAp ' 

which finishes the proof of the sublemma. □ 



Sublemma 3.5. Let j ^ 1 he an integer, let pj and aj be the real and imaginary parts 
o/log(l — yZ/Sij-ie*^), where is a beta random variable with /3{l,j — 1) distribution 

and 9 is independent of uniform on [0,27r]. Then, for s,t G M, 

if s'^ + t^ ^ Sf, and 

Q 

|]E[g«(*Pj+'5o-j)j| ^ 

if s'^ + t^ ^ 8j2 and i ^ 2. 



+ 12 



Proof. For t G M and s G C with real part strictly between —1 and 1, 

Tgr r(j)r(j + it) 

^ ^ r(j + ^)r(j + ^) ^''^ 

(see [3]). Now, if t is fixed, the function 

is holomorphic, since the imaginary part is uniformly bounded (by 7r/2), which implies 
that (fTO|) holds for alH G M, s G C, and in particular for all s, t G M. Moreover, 

T{k)T{k + it) ^ 

r (fc + ^) r (A; + ^) fe"^ ' 

since r(fc + z)/T{k) is equivalent to for all 2; G C. Hence, by using the equation 
T{z + 1) = zT{z), one deduces: 

If + ^2 ^ 8j2, then |Q(fc,s,t)| ^ 1 for all k^ j and |Q(A;,s,t)| ^ e-(^'+*')/i°^' for all 
k ^ 3j. Hence 



00 

k=3j k=3j 
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Now let US assume + ^ 8j^. One has: 

^ J r(i + ^)r(i + ^) llQ(fc,s,t) 

where all the factors q^^}^ have absolute value bounded by one. By considering the case 
where j = 1, one deduces 



r(i)r(i + zt) 



r(i + ^)r(i + ^) 



^ 1, 



and then, for j ^ 2, 



\Qil,s,t)\ Vi^Tt^ 

□ 

Sublemma 3.6. For N A and 9 G [0,27r), the distribution of \og{Zx{0)) under Haar 
measure on U (N) has a density with respect to Lebesgue measure on C, which is continuous 
and bounded by Cq/ log{N) , where Cq > is a universal constant. 

Proof. By the results in [3] and the previous sublemma, one checks that the characteristic 
function $ of \og{Zx{9)) G C ~ is given by 

N 

If + ^2 ^ 32A^, one has + ^ 128 ^ 8f for j G {2, 3, 4}. Hence, 



s 



+^2)3/2- 



If s^ + t^ ^ 32N, then s^ + t^ ^ 8f for all j ^ Hence, 



ms,t)\<: H E[e^(*''^+-^)]^exp(-(.2 + t2) ^ _L j 

2v^s£KAf \ 2VlVs£KAf / 



Since e^/-' > one deduces 



N+1 



3VN 



(s2+t2)/30 



^ I _J_ I ^ g-log{Af/9)(s2+t2)/60 
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Now, for ^ 10, 



\Jo Jz2N U ' J 

GOtt ^ , 1/2 10000 

^ — 7T7TT + 10247r 32Ar ^ 

log(A^/9) ^ ' log iV ' 

and for N e {4,5,6,7,8,9}, 

f 512 f 

\^{s,t)\dsdt ^ / --^——^-^l{s2+t2^32N}dsdt+ / l{^2+i2^32Ar}ds(it 

JK2 [S +t ) ' J^2 



r32N roc 

nil du+ I ^-3/2 ^"^ 



512 

J32N " 



^ 327rA^ + 10247r(32iV)-^/2 ^ 2887r + 102471(128)-^/^ ^ 10000^ 

log 9 

By applying Fourier inversion, we obtain Sublemma 13.61 □ 



Let us now go back to the proof of Lemma 13.31 For any X G U{N) with eigenvalues 



Qj^g j-^g^g^ pg^gg -y^^j^grg ^ g«fj for all j G {1, . . . , A^}, and modulo tt. 



X := Im(log(Zx(e)) = J] Im(log(l - e'^''^-^))) 

= i ^ (0^. _ 0) + ^ Im(log(e-^(^^-^)/2 _ g^(e,-e)/2)) 

= ^Im(logdet(X))-^+ J] Im(log(-2*sin(^,-0)/2))) 



2 Im(log det(X)) - ^ 



where JT" denotes the version of Im(log det(X)) lying on the interval (— vr, vr]. Hence, for any 
e G (0, vr), \ J'\ ^ e if and only if X is on an interval of the form ^ 2kTT+e-N {e+n 

for some G Z. Now, for some A > Q chosen later in function of 5, let $ be a continuous 
function from C to [0, 1] such that ^{z) = 1 if | Re z—xq \ ^ 6^/[ogN and | Im 2;| ^ A^/logN, 
and such that ^{z) = for | Re^ - xo| ^ 25^/\ogN or | Im^l ^ 2Ay/\ogN. For e G (0, vr). 
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and under the Haar measure P;7(Ar) on U{N), 



= — ^E[/(Ar) $ (lOgjZx (^) ) ) 1 r 2k7.-e-N(0+n) 2k7: + ,-N{e + 7:) j 

feez 

= — / dx dy D{x + iy)^{x + iy) 

/oo />l/2 
dx / dML'(x + i[A;7r-iV(^ + 7r))/2 + 'ue])$(a; + i[A;7r-7V(^ + 7r))/2 + «€]), 
-oo J-1/2 



TT 

feez 



-1/2 

where D denotes the density of the law of log(Zx(^)), with respect to the Lebesgue measure. 
Now, 

D{x + i[kn - N{9 + 7r))/2 + ue])^x + i[kTT - N{9 + 7r))/2 + ue]) 
is uniformly bounded by the overall maximum of D and vanishes as soon as |a; — XqI ^ 
25 Vlog N or |A;|7r ^ iV(|6'| + 7r)/2 + 7r/2 + 2A^/^ogN. Since D and $ are continuous 
functions, one can apply dominated convergence and deduce that 

^EuiN) [$(log(Z^(^^)))IL{|^|^,}] 

converges to 

/oo 
D{x + i[kn - N{0 + 7r))/2])$(x + i[kn - N{0 + n)) /2])dx 
... .. -oo 



when e goes to zero. On the other hand, if the matrix X follows '^sv{n) and if T is an 
independent uniform variable on (— tt, tt], then Xe^'^l^ follows P[/(Ar) and its determinant 
is e*"^. One deduces: 

-^V(^N) [$(log(Z^(^)))l{|^|^,}] = -^SV{N) [$(log(Z^,,T/.v(^)))]l{|r|^e}] 

= 1 y Esvi,N) mog{Z^^^i^{Q)))\ dt 

1/2 

^SU{N) Ml0g{Zxe^i,e/N{e)))] dv 

1/2 

Now, the function X i— )> $(log(Zx(^))) is continuous from U{N) to [0, 1], since $ is contin- 
uous with compact support and X i— )■ log(Zx(6'))) has discontinuities only at points where 
its real part goes to — oo. One can then apply dominated convergence and obtain: 

^Eu^r,) [$(log(ZxW))l{|^K.}] ^^^su(N) miog{Zx{em. 

By comparing to the convergence obtained just above, one deduces 

/oo 
D{x+i[k7r-N{0+7r))/2])^x+i[k7r-N{0+7r))/2])dx. 
... .. -oo 
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Since D{z) ^ Co/logN and 



■{|a;-xo|!j<5yi5i7V,|j/| 

for all x,y eM, one deduces 



SUiN) 



ndLCo 
logN 



where d = A6^y\og N is the length of the interval [xq — 26^/\ogN, xo + 26^/\ogN] and L is 
the number of integers k such that \k7T — N{6 + 7r))/2| ^ 2A^/[ogN. Now, it is easy to 
check that L ^ 1 + ^A^/ME and then 

77 ' 



P 



SU(N) 



log|Zx(^^)|-Xo| ^5v/b^,|ImlogZx(^^)| ^ Av^b^ 



^ 16 Cq Ab + 



47r5Co 
VlogiV' 



Using Lemma [3 ■2[ one obtains 



P 



SV{N) 



log|Zx(^^)|-xo| ^(5v/b^ 



^ 16 Cn Ab H + c^t 



^J\^gN 



Let us now choose A := 1 + 51og(l/5). One gets 



and then 



Therefore, 



, A A 
2 V 2 



2 ' ov , .J 2 

v4^\ ^ 5v^log(l/5) 



^ log(l/5). 



P 



SU{N) 



ATTbCn 



\og\Zx{e)\ - xo\ ^S^logN ^ 16Co6 + 80Co6log{l/6) + -=^ + c[6. 



VhgW 

Since 6 < 1/2, one has 6 ^ (51og(l/5)/log(2), which implies Lemma [3. 3[ for 



C 



16 Co 
log 2 



80Co + 



47r Cn 



+ 



(log 2)3/2 log 2 



□ 



3.5. Behaviour of the oscillation in short intervals of the log-characteristic poly- 
nomial. 

Lemma 3.7. There exists C2 > such that for /i G M and A ^ and uniformly in 
^ M ^ 2 V 



P 



SC/{Af) 



P, 



S(7(Af) 



RelogZx + ^) - RelogZxiO) ^ ^ AVlogM^ ^ ^ 
^ ImlogZx + ^) - ImlogZx(^^) ^ ^ AVlogM^ ^ 



C2_ 

A^' 
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Proof. By symmetry of the problem, we can assume fi > 0. Setting 

Re := Re log Zx + ^) - Re log Zx{e) 
for fixed (or the same expression with the imaginary part), we get: 



(r dO \ 1 / / f dO \ 



1 


log 


;M 


1 




AHog 


;M 


1 




log 




1 





2" ,2C/0\ 
\S.SUiN)[Rg) — 

EuiN){Rl) (by m) 



A^ log M 

Now, under U{N), the canonical matrix X is almost surely unitary: let 6i,. . . ,6n be 
its eigenangles in [0,27r). For j G {!,..., A^} and t G [0, 27r)\{6'j}, we can expand the 
logarithm: 

log(l-e^(^^-*)) = -5^^^ 
as a semi- convergent series. Hence, for t such that Zx{t) ^ 0, 



A; ^ k 

j=l k^l k^l 



Thus: 



RelogZ.(t) = -i le--tr (X^) + le^tr (X^)) = -'-J^ ^e'-tr (X^) 

Im log (0 = -i f E r - E r'" (■^") ) = E r (■^■') ■ 

Here, the series in A; G Z* are semi-convergent: more precisely, setting for K ^ 1, 

ke:Z*,\k\^K ' ' 

and 

St := RelogZx(t), 

S'}^'' tends almost surely to St when K goes to infinity. 

Moreover, one has the following classical result (^): for all p,q ^ "Z, 

EuiN) (tr (Xf )t7IX^) = \p\ A X. (11) 
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Hence, for i^T, L ^ 1, t, n G M, 

/ \ / 1 p—i{pt+qu) 

\ p,qeZ*,\p\^K,\q\iiL 
1 -i{pt+qu) 

p,q&*,\p\<K,\q\^L 
X ^-iipt+qu) 

= - Yl ^^1{P=-,} kl A AT (from (in 

= 2 ^ -^cos(/c(M-t)). 

One deduces that 

E^(^) ((5f ) - S['^f) = E^(^) ((5l^V) + E^(iv) ((^i^y) - 2 E^(^) (^i^^Sf 



1 A AT , ^ 

9 ,2 ^'^^ ^^^^ " l{i^AL<fc<i^VL} 



2 ^ fc2 



A; A AT 



2 ^ A;2 

fcSKAL 



which tends to zero when K AL goes to infinity. Hence, Sj.^^ converges in when K goes 
to infinity, and the hmit is necessarily 5*^. Therefore, 

^ /r, o N , 1 A;AA^ , ,, l>r^fcAA^ , ,, 

^u{N){StS^) = hrn^- — p— cos(A;(M-t)) = -2^— p— cos(A;(M-t)). 

The same computation with St := ImlogZx(t) gives exactly the same equality: 

A; A 



It is therefore enough to achieve the computations only with St- Using this last formula. 



we can write, with a = -^i 
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^U(N){Ro) = ^U{N){{Sa — Sq)'^) = 2Eu[S'^ — SaSo) 



We can then develop E[/(7v)(-Ro)- 



^u{N) [Ro] = 2^ — - 2^ cos 



k>l 



k^ 



kfi 



But we also have 
kAN 



k^l k=l k>N ^ \ / / 

Moreover we have the following result ([7j, p. 37), uniformly on ^ G [— 7r,7r]: 



2^— p— cos(A;6') = - log 2 sm I - 

k>l ^ 



TT 



+ Ci (A^ 1^1) + cos (A^^) --N\e\+Ne Si (A^^) 



+ 



(12) 



where: 



Si iz] 



smx , vr cosz 

ax = h 

X 2 z 



cosx 



dx 



Ci (z) 



cosx 



X 



dx = ■J — log z 



x^ 

cos X — 1 

X 



dx 



Recall also that 



smx vr 

dx = —. 

X 2 



Let us denote / (/i) := log/i + |/i — cos/i — Ci (/i) — fiSi (/x). We have, for A^ going to 
infinity: 



^u{N){Rl) = logAr + l + 7 + /(/i)-log/i + log 



2 sin 



0(1/A^) (with 



log AT + 1 + 7 + /(/.) - log/. + log (2^ (^1 + O ((^)'^^^ + 
l + 7 + /(/.) + 0((^)^) 



0(1/Ar) 



0(1/Ar). 



Let us now study the behavior of the function /. 
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fin) = log/x-cos/x + /x -Si(/x)) -Ci(/x) 



/ cos fjL 



, , „ , COS X \ f , COS X — 1 . 

log IJ, — COS + /X ( / 2~ ) " ( 7 ~ + / 1 



COS a; , cos x - 1 , 
= —7 — n — 5— aa; + / ax 



IJ, Jo 



One has: 



= -7 + 0(l) + 0(l + log(//Vl)) = o(log(^V2)), 



which imphes 



uiN) {Rl) = O (log (A V 2) ) = O (log M) . (13) 

□ 



3.6. Control in probability of the mean oscillation of the log-characteristic poly- 
nomials. 

Lemma 3.8. For a certain n e N, let us consider an i.i.d. sequence {Uj)i<^j<^n of random 
matrices following the Haar measure on SU (N) . Let us set: 

y'iiogjv 

For 5 e (0, 1/2), let us consider the random set: 

Ss ■= |j({^e[0,27r] / |L,(^)|^ri}u|J{^e[0,27r] / |L,.(^)-L,(^)|^(5}) 



Then, there exists C3 > 0, depending only on n, such that for all N A: 



E(A2.(<f5)) ^ C35l0g(l/5), 

where A27r denotes the normalised Lebesgue measure on [0,27r]. 
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Proof. Using Markov's inequality, one gets: 



Now, 
and 



27r 



d9. 





1=1 ''^ l^i^j^n-'^ 

— P5C/(JV) \\\og\Zx{e)\\^5-'^-\ogN 
+ n{n-l) sup ¥su(N)\\\og\Zxm-A^5\\^ogN 

ee[0,27r],xGlR \ V ^ 

^ ncie"^(V^'^) + n{n - l)C{5/V2) \og{V2/5) 



e ^V^" 2 V2\V2'' 2 / ^ e"— = 0((51og(l/5)) 

51og(x/2/5) ^ 51og(x/F^/5) = y log(l/5). 



which gives Lemma [3.81 □ 



3.7. Control in expectation of the oscillation of the log-characteristic polyno- 
mials on a small period. In the sequel, we consider the dimension ^ 4, an integer K 
such that 2 K ^ ^/2, defined as a function of N which is equivalent to A^/(logiV)^/^^ 
when N goes to infinity. We denote 

M := N/K ^ 2, 

which is equivalent to (log A^)^/^^, and we also define a parameter 5 G (0, 1/4) as a function 
of A^, equivalent to (log A^)"^/'^^ when A^ goes to infinity. For ^ [0,27r], we denote, for 

27ck 27ckM 
"k '■= 00 H — — = 00 + 



and for ^ ^ - 1, 



K ' N 



A := 0ls■J-^ — 01- = — = . 



The angle is chosen in such a way that the following technical condition is satisfied: 

K-l 

^SUiN) 

d0 



k=0 

^ KE 



SU{N) 



(|lmlogZx(^fc + (1 - v^)A) - ImlogZx(4 + v^A) 

ImlogZx(^ + (1 - V5)A) - ImlogZx(^ + v^A) ) . 
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This choice is always possible: indeed, if the converse (strict) inequality were true for all 
6^0, then one would get a contradiction by integrating with respect to 6*0 G [Q,2t[/K). We 
then define the interval J := [6*0, Oq + 27r) = [6*0, 6k)- Note that all the objects introduced 
here can be defined only as a function of N . Moreover, by applying Lemma [3?71 to 6* + v^A 
and /i = A^(l — 2a/5)A ^ 27rM, we deduce that the assumption made on 9q implies: 



K-l 



J2 ^suiN) (|lmlogZx(^fc + (1 - v^)A) - ImlogZx(4 + v^A)|) = O (^K^hgM^ . 

k=0 

(14) 

We can then introduce the 2-oscillation of the real and imaginary parts of the log- 
characteristic polynomial: 

Definition 3.9. For 9 E J and fi E [0, 27rM], and for the canonical matrix X E U{N), the 
2-oscillations of Re log Zx and Im log Zx are defined by 



A^/e 



VMM) 
1 



Re log + ^) - Re log Zx (9) 

lm\ogZx{9 + ^^-lmlogZxi9) 



In case of several matrices {Xj)i<^j^n, we denote the corresponding 2-oscillations by 
A^Ri'^ and A^4^l 

In the sequel, we need to introduce several random sets. The most important ones can 
be informally described as follows: 

(1) A set ^ of indices k such that the average of the 2-oscillations A^Rq and A^/g of 
the log characteristic polynomials for 9 E [9k,9k+i\ and E [0,27rM] is sufficiently 
small. 

(2) For k E ,JY\^ a subset of [6*^, 9k+{\ for which the average of the 2-oscillations with 
respect to /i G [0, 27rM] is small enough. 

(3) A subset ^ of jV\ of "good" indices, such that there exists 9 E \9k,9k^ v^A], 
both in and (f/. This last set, introduced in Lemma [3. 8 [ corresponds to the fact 
that the logarithms of the absolute values of the characterize polynomials are not 
too large and not too close from each other: from this last condition, we can define 
the " carrier wave" . 

(4) For k E and for some 91 E [9^, 9k + v^A] n % n <ff^, a subset of [0, 2'kM] 
such that the 2-oscillations A^Rg*^ and A^Ig* are sufficiently small. This condition 
will ensure that the carrier wave index corresponding to 9 = 91 + fi/N does not 
depend on G 

(5) From this property, we deduce that, for each pair of consecutive gaps between zeros 
of the carrier wave, which are sufficiently large to contain an angle of the form 
91 + fi/N for k E ^2 and fi E Wk ("roomy gaps"), one can find, with the notation 
of the introduction, a sign change of i^e*^^/^F/v(e~*^), and then a zero of F^. 



SUMS OF CHARACTERISTIC POLYNOMIALS OF UNITARY MATRICES 25 

All these sets will be precisely defined in the sequel of the paper, in a way such that their 
measure is "large" with "high" probability. The corresponding estimates will then be used 
to prove our main result. 

Lemma 3.10. Let be the n-fold product of the Haar measure on SU{N), E^'^^jy^ the 

corresponding expectation, and {Xj)i^j^n the canonical sequence of n matrices in SU{N). 
Then: 

(1) There exists a random set ^ C |0, -ft' — 1] such that E^'^^^^(|o/(^|) ^ (1 — 6)K and 

p(n) 



^^"^(^)-a.s.,V(j,A;)e|l,n]x^, 

^^"'^ 27rM27r " \6K 



'k+l 



and 



1 io " 27rM27r \5K ) 



(2) Pl?f)(^)-a.s., \/k e ^i, 3% C [9k,0k+i) such that X2n {%) ^ (1 - S)/K and, 

""a R^^^^ - O and r" A J^^)^ - O ( 1-] (15) 

i^ere, i/ie implied constant in the O(-) symbols depends only on n. 

Proof. By f lT3|l and the similar estimate for the imaginary part, we have uniformly (with 
a universal implied constant), 

Ef;(;v)((A^i?o)') +Ef/(^)((A^/o)') = 0(1). 

The Cauchy-Schwarz inequality ensures that 

Ec/(^)(A^i?o) + Ef;(;v)(A^/o) = 0(1), 

i.e. 



dd 

E5^(^)(A^/?, + A^/,) — = 0(1), 



J 27r 



which implies 

II W.(A.i^« + A„/.)5^-^0(l). 

Splitting the interval J into K equal pieces and applying this estimate to n independent 
matrices (Xj)i^j^„ following the Haar measure on SU{N), one gets 

^EE-f^,4f"r'(AX''.A„tf')^|)^0(l,. (16, 

j=l k=0 \JOk -JO / 

Applying Markov inequality, we deduce that there exists a universal constant k > 0, such 
that 
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Eg)(^) (^card I (j, fc) e II, nl X fO, K-l}/ j^^ 
and 

Eg}(^) (^card I (j, fc) e |1, n] X |0, - 11/ ^ 



''^^ 2^2^ ^ K5 



and 

I ^ov^ J ^ iTi ,K,n iTn _ m / / / /> 

7o ^ 27rM 27r ^ ^5 

We thus set 

L 7o ' 27rM27r ^ 



and we get: 



4™(i=^i)^(i-<^)^- 



Now, for k e let us set: 

Applying again Markov inequality, we get that P^^^^^-a.s.: 

A2. {%) ^ (1 - S)/K. 

We now define good indices. 
Definition 3.11 (Good indices). An index k E {0, K — Ij is said to be good if 

(1) ke^K, 

(2) ^in%n[9k,9k + V5A)^$ 

We denote by ^ the set of good indices : 

^2 := {a: e ^ / 4^ n n [9^, 9^ + V6A) + 0} 

An index is said to be bad if it is not good. 
We then get the following result: 
Lemma 3.12. With the notation above, the set of good indices satisfies: 

^suiN)i\-^2\) =k{i-0 (v^log(l/5))) , 
where the implied constant in the O(-) symbol depends only on n. 
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Proof. If k e ^2"^ either k G or A; e ^ and n%n [9k, Ok + v^A) 
J/^"" = U where: 



I.e. 



^i:={ke^i 



By (HZD, we have E^g^^^d^^]) ^ 

For all /c G J?^ we have ^5 D ^fc H [^fc, + v^A), i.e. ^ D Ufce^iT^^ ^ t^^'' + ^^)' 
where the union is disjoint, and thus, \2-k ((^5) ^ niinfc A27r ( '^k H [^^a,., ^^a,. + V^A) ) . 



P 



(n) 



SU{N) 
A2 



•a.s., we have: 



^(^fcn[4,^fc + v^A)) ^ A2.(?^fc) + A2. ([^fe,^fc + v^A)) -A2.([^fc,^fe + A)) 



Now, since 5 < 1/4, we obtain P^^j^^^-a.s.: 



A2. (^A.n[^fe,^, + v^A)) 



(20) 



This implies that P^j-^-j-a.s.: 



Now, by Lemma (13. Sp . E^'^^.^^ 



O (Kv^log(l/5) ) and then: 



E 



(n) 



^ 5ir + f/s:v^log(l/(5) 



□ 



3.8. Speed of the good oscillation of the log-characteristic polynomials. 

Lemma 3.13. With the notation above, and P^j-^^-a.s., VA; G there exists a random 
set ^k C [0, 2ttM], and 61 G n%n [Ok, 0k + VSA), such that 

A2.M m) = 1-0 (r2(logiV)-V4(iogM)i/2) , (21) 
where Xm is 1/2ttM times the Lebesgue measure, and for all j G [l,?^], /U G Wk, 



(log M) 1/2; ^ \{log My/^ 

Again, the implied constant in the O(-) symbol depends only on n. 



m 
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Proof. Let ke^2 and 61 E n n [Ok, 9k + VSA). We set: 

n 

^ := n {A.i?|^ ^ e, A.4f ^ e} 

where 

_ (log AT) 
^ (log M)i/2- 

Applying Markov inequality, we get: 



e i^jX" Vio 27rM Jo 27iM J e ^ ' 



by (fTSjl which gives the announced result. 



□ 



3.9. The number of sign changes. 

Let us go back to Theorem 11.11 We need to estimate the number of zeros of on the 
unit circle, or equivalently, the number of values oi 9 E J such that the following quantity 
vanishes: 



i=i j=i 

Using the fact that U^j E SU{N), one checks that e^^^^"^ Zu^ _.{9) is real, and then the 
number of zeros of F/v on the unit circle is bounded from below by the number of sign 
changes, when 9 increases from 9q to 9q + 27r, of the real quantity given by the right-hand 
side of fl22l) . Now, the order of magnitude of log \ Zuj^ {9)\ is ^/\ogN and more precisely. 
Lemma 13.81 informally means that for most values of 9, the values of \og\Zu^ - {9)\ for 
^ ^ j ^ n are pairwise separated by an interval of length of order y/\ogN . Hence, one of 
the terms in the sum at the right-hand side of (122|) should dominate all the others. If j is 
the corresponding index, one can expect that the sign changes of (j22l) can, at least locally, 
be related to the corresponding sign changes of i^e'^^^^'^Zuj^ .(9), which are associated to 
the zeros of the characteristic polynomial Zu^.. This should give a lower bound on the 
number of sign changes of f l22l) . 

This informal discussion motivates the following definition. 

Definition 3.14. With the notation of the previous subsections, for all k E J^2-, we define 
the carrier wave index by: 

ik ■■= Arg max {Re log Zx^ {91) ] , 
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where 6'^ is the random angle introduced in Lemma 13.131 Moreover, we consider the 
following interval: 



Jk 



As ei e 4^ we have Vj ^ jk, RelogZxM'- 
we deduce that Vj 7^ j^, V/i G : 



^ RelogZx,._(^^) - ^VWN. From (|2l 



RelogZ^^, (ei + ^) ^ RelogZx,^ {ei + ^) - ^(logiV)^/^ + O ((logiV)^^) (23) 
Now, since 

l/5 = 0((logiV)ViO)^ 
with a universal implied constant, we then get, for a universal c > 0, 



^ exp (-2c(log Ar)°-^ + O ((logA^)^/^)) ^ exp (-c(log Ar)°-^) 



for large enough, depending only on n. This implies: 



miuj \bj\ 



ii 



exp (-c(logAr)°-^) 



N 



for A^ ^ A'o, where A'o depends only on n,bi, . . . ,bn- Hence, for k G E '3^k and 

= 91 + fi/N, the quantity 

which is P^^^^^-j-a.s. real, has the same sign as its term of index jk- 

Theorem II. II is proven if we show that the expectation of number of sign changes of G{9) 
for 6 E J, under P^.^s, is bounded from below by A^ — o{N). Hence, it is sufficient to get: 



'^SU{N) 



( ^N-o{N), 



where Sk is the number of sign changes of hj^i^ e'^^/'^ Zx^^ (0), for 9 e Jkn{9l + fi e ^}. 

Now, for k G let ak,i ^ ak,2 ^ ■ ■ ■ ^ (^k,^^ be the eigenangles, counted with 
multiplicity, of Xj^, in the interval Jk- The sign of bji^e^^^^'^Zxj^ alternates between 
the different intervals (aA:,i, «a:,2), («fc,2, «fc,3), • • • , («fc,i/fe-i, ctA:,!/^)- Hence, for each pair of 
consecutive intervals containing an angle 9 = 91 + -^jfi G we get a contribution of at 
least 1 for the quantity Sk- 
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Every element of Jk can be written as 6^^ + for 

^ /i ^ (1 - v^)iVA ^ iVA = 27rM. 
The Lebesgue measure of the elements of Jk for which /i ^ is then bounded by 



2'kM 



N 



where A denotes the standard Lebesgue measure. Hence, if an interval {ak,u, has a 

length strictly greater than this bound, it necessarily contains some 6* = + for which 
/i G For some c' > depending only on n, this condition is implied by 

- ak,. > c'^5-2(logiV)-i/^(logM)V2. 

We will say that (a^^jy, <yk,u+i) is a roomy gap if this inequality is satisfied, and a narrow 
gap if 

ak,u+i - ^ c'^6-\\ogN)-'/\logM f'. 

By the previous discussion, Sk is at least the number of pairs of consecutive roomy gaps 
among the intervals {ak,i, ak,2), {ak,2, "fc.s), • • • , {c(k,Uk-i, C(k,uk)- If there is no narrow gap, 
the number of such pairs is (z/^ — 2)_|_ ^ t'fc ~2. Moreover, if among the intervals, we replace 
a roomy gap by a narrow gap, this removes at most two pairs of consecutive roomy gaps. 
Hence, we deduce, for all k G that 

5fc ^ z/fc - 2 - 2tpk, 

where Uk is the number of zeros of Zx in the interval Jk and ipk the number of narrow 
gaps among these zeros. Hence, we get the lower bound: 



E 



(n) 

SU{N) 



^ E' 



in) 

SU{N) 



( 5^z/fc-2ir-2^j 



where is the total number of narrow gaps among the zeros in [0, 2tt) of all the functions 



Now, P 
J^k 



in) 



-a.s., for all k G we have: 



' SU(N) 

{ee [ei,ei + {i-V~6)A],z,,{e) = o} 
> {de k + v^A, dk + {i- 75) aI , Z.j^{d) = 0} 



= + 1 (im log Zx, (^fc + (1 - v^) A) - Im log Zx,^ (Ok + v^A) 

N 1 " I 

^ _(1 _ 2v^) - - ^ I (lmlogZx,{ek + (1 - V5)A) - ImlogZx^iOk + VSA) 

the second equality coming from Proposition 12.21 
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Adding this inequality for all k G taking the expectation and using (iMj) yields the 
estimates: 



^SUiN) 



n K-l 
j=l k=0 



(n) 

SU{N) 



N 



ImlogZxM + (1 - V5)A) - ImlogZv.(^fc + v^A) 



^ —(1 - 2V6)K{1 - 0{V6\og{l/6)) + 0{K^/\ogM) 
K 

>iV(l-0(V?log(l/^))+of^^^ 



Moreover, 

It remains to estimate 



2K = 0{N/M). 



(24) 
(25) 



4^(7V)[2V'] = 2n¥.su{N)[x] = 2n¥.uiN)[xl 



where x denotes the number of narrow gaps between the eigenvalues of the canonical 
unitary matrix X. The replacement of SU{N) by U{N) is possible since the notion of 
narrow gap is invariant by rotation of the eigenvalues. 

Now, the last expectation can be estimated by the following result: 

Lemma 3.15. For N ^ 1 and e > 0, let U be a uniform matrix on U{N) and let Xe be 
the number of pairs of eigenvalues of U whose argument differ by at most e/N . Then, 
nXe]=0{Ne^). 

Proof. For ^^i, ^2 ^ I^, the two-point correlation density of the eigenvalues of U at e*^^ and 
e*^2, with respect to the uniform probability measure on the unitary group, is given by 



p(e 



Now, 
and then 



&m[N{d2 - di)/2] 
Arsin[(^2-^i)/2] 

A^|sin[(e2 -^i)/2]| ^ N\e2-ei\/2 

sin[Ar(g2 - gi)/2] " ' 
Arsin[(^2-^i)/2] 



smx 



X 



for X = N{92 — di)/2. Now, for all x G M, | sinx| ^ sin |x| ^ |x| — |a;|^/6, which implies 



smx 



X 



^ 1 



X 

"6 



2\ 2 



> 1 - 



X 
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1 - 



/ sinx 
\ X 



6 



Integrating the correlation function for 9i G [0,27r) and 0' := 62 — 6i E [—6/N,e/N] gives: 

o c J-e/N 



2tt J_^/n 271 6 



□ 



From this result, applied for 



e = c'M6-\logN)-^/\logM) 



vl/2 



we get the estimate: 



E' 



(n) 

'SU{N) 



0{Ne^) = O {NM^6-%logN)-^/\logMf/^) . 



(26) 



The estimates ([24]), ([25]) and ([26]) imply: 



'^SU{N) 



^ N 



1-0 { v^log(l/5) + ^^^^ + M3r6(logiV)-3/4(iogM)3/2 



M 



From the values taken for 6 and M, we get: 

V6\og{l/6) = O ((logiV)-=^/^MoglogAr) 



yiogM 
M 



O (^VloglogArQogiV)-^/''^ 



and 



M=^r6(l0gAr)-=^/^(l0gM)3/2 = O ((logAr)9/64(log]V)18/32(logiV)-3/4(loglog]V)3/2) 

= O ((logAr)-=^/^^(loglogAr)3/2) . 

Finally, we get 

An) 



E 



SU{N) 



J2^^=N{l-0{i\ogNr/^^)) 



which completes the proof of Theorem ll.li 
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